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SINOPSIS 


Thesis entitled ’Some Studies in Dilepton Production in 
Hadron-Hadron Collision' by Mohammad Noman in partial fulfil- 
ment ;of the requirements of the Ph.D. degree of the Department 
of Physics, Indian Institute of Technology, Kanpur. 


^Quantum chromodynamics (QCD) seems to be a suitable 

theory for strong interactions. It is a non-abelian gauge 

theory of interacting coloured quarks and gluons. The quarks 

come in three colours and the gauge fields called gluons form 

an octet under SU(3), the underlying symmetry being STJ(5) . 

Colour 

-It is the asymptotic freedom that maKes QCD interesting and 
applicable to high energy processes. The effective QCD coup- 
ling constant goes to zero logarithmically as the momentum 
transfer squared goes to infinity, i.e. at very high momentum 
transfer, the theory becomes effectively a free field theory 
and hence perturbation theory can be applied in such processes. 
Many calculations have been done involving hadrons in the 
asymptotic region of QCD and there seemb to be fairly good 
agreement between theoretical perturbative calculations and 
the experimental data, such as - 

(i) scaling violation observed in deep inelastic lepton- 
hadron scattering, 

(ii) the ratio R = the e^e~ annihila- 

a(eV- Y* .|iV) 

tion, 
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(iii) hadron-hadron hard scattering producing particles with 
large transverse momentum 

(iv) jets of hadrons observed in 6*^6” aruiihilation and other 
hard scattering processes, and 

(v) massive dilepton production in hadronic collisions. 


.In this thesis we investigate some more aspects of 
dilepton production in the QCD framework. We study the angular 
distribution of the dilepton in their rest frame, the z axis 
being along the beam direction and x axis in the plane of the 
colliding hadron's momenta, i.e. the Gottfried- Jackson frame.. 


The QCl subprocesses considered are qq -^gy " ^ _ and qg'^qy 


* 




The angular distribution is of the form 1 + A cos 6 + B sin 29 
cos 4'+ C sin^e cos 2({>. A, B and C = are functions of 

dimuon mass m, s and q^, (the transverse momentum of the dimuon) 
,We neglect the primordial transverse momentum of the partons, 
and restrict ourselves to large q,j^, in order to avoid mass 


singularity and infrared divergence.- ¥e compute A and B as a 

% - + 

function of r^ = ^ for pp, pp, up and n p collisions at 
ifg = 27.4 Gev. ' A and B show striking variations with r^. ¥e 
find that for low values of r^ '^0.05^ the angular distribution 
is of the form 1 + cos^ ® (A - 1, B = C O) identical to the 
Drell-Yan prediction. However, for higher values of r^> 0.05, 
there is a strong dependence on both 9 and <!> , which shows a 
striking departure from the Drell-Yan prediction. For example, 
for x = m^/s = 0.1 and r^ = 0.3, A = 0.18 and B = 0.12 for pp 
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collision. Hence it is the high transverse momentuni region where 
the QCD shows its effects without any ambiguity. 


Purther we study the asymmetry 

da(H3_(+)H2(+) ^ - do( H^( +) -) 


^ do(H3_(+)H2{ +) -^ + do( H^( +) H2( -) 

associated with dilepton production in longitudinally polarized 
hadron-hadron collisions H-, +Hp'*' p'*'p"’x. ¥e study the variation 

'll ^ 

of Aj^ with ±n longitudinally polarized p-p and p-p 

collisions using four different models of spin dependent quark 
and gluon distribution functions. The basic QCD subprocesses 

For p-p collision 


— ■X' -)£. 

considered are qq -^-gy ^ _ and qg->qy ^ _■ 

p p, p 

the asymmetry is quite large (greater than 15?^) and shows striking 


variation with r^ and is strongly dependent on the model of spin- 
dependent structure functions. For pp collisions the asymmetry 
is not large ('^5 to 7fo) and it is not strongly model dependent 
because in all the models, we use, the sea and the glue are not . 
strongly polarized. -The study of the asyonmetry parameter Aj^ as 
a function of r^, serves two purposes:* (i) it distinguishes the 
QCD subprocesses giving rise to dimuon where Aj^ will be non-zero 
from the quark-meson scattering- subprocesses giving rise to. 
dimuon for which Aj^ will be zero and (ii) it serves as a good 
test for the spin dependent quark and gluon distributions in 
polarized nucleons. 

* Finally we derive relations between dilepton differential 
cross-sections, based on 1st order QCD subprocesses, in meson- 
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nucleon collisions: Tt— N, k— N, k°N, k'^N, {N = n,p) collisions, 
assuming the sea to he SU( 2) singlet. These relations lead to 
tests of the parametrizat ions of (i) nonstrange x^^alence quark 
distribution in mesons and (ii) gluon distribution in n and k 
mesons? The knowledge of the sea quarks either in mesons or 
nucleons and gluon distribution in nucleons is not required. 



Chapter 1 
INTRODUCTION 


Quantum chromodynamics (QCD) seems to he a suitable candi- 
date for strong interactions. It is a nonahelian gauge theory of 
interacting coloured quarks and gluons. The quarks are supposed 
to come in three coloxars and the intermediating gauge fields 

called gluons form an octet under SU(3), the underlying symmetry 

1 2 

being SU( The QCD Lagrangian ’ is 

- I - ’v + Sf®"” 1,-iS r ♦ 

^ t' 

( 1 . 1 ) 

where are the 8 gauge fields, ^ is the quark field, g is the 

T_ _ 

coupling constant, f®' are the structure constants of SU( 3) and 
are the generators of SU(3) in the triplet representation. 

In contrast with quantum electrodynamics, the gauge fields 

of QCD carry colour and hence interact among themselves. The 

Reynman rules^ for QCD are summarized in fig. 1.1. Although 

there are many hurdles such as the problem of quark and gluon 

confinement etc. before QCD, yet it becomes useful in calculating 

processes involving hadrons at very high energies because of 

1 2 

'asymptotic freedom'. ' 


By 'asymptotic freedom' we mean that the 
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2 

coupling constant decreases with increasing Q , the momentum 
transfer squared. This remarkable result of asymptotic free- 
dom arises in the solution of the renormalization group 
equation for n point Green function of non-abelian gauge 
theories. One defines a running coupling constant g(g,t) 
(where t = In Q^/Q^) by 


1 p(g); g(g,o) = g 

p(g) is the so-called Callan-Symmaxizik function. Por QCD 

.3 


( 1 - 2 ) 


P(g) 


4n 


b . (33 - 2f) 


(1.3) 


and fortunately p(g) is negative if the number of quark flavours 
f 5. 16 and hence g(g,t) + 0 logarithmically as Q -► ® - Eqs. 
(1.2) and (1-3) immediately give 


r(t) = 




1 - li i"(o) 


(1.4) 


In deep inelastic hadronic processes, it is assumed, 
because of aymptotic freedom, that the hadron constitutes of 
almost free quarks and gluons and these quasi free partons 
interact incoherently with the high Q probe (virtual photon 
etc.) and hence one needs quark and gluon distribution func- 
tions for processes like hadron ^ parton and quark and gluon 
fragmentation functions for processes like parton-*- hadron. 


In these processes two time scales'^ are involved. One is a 
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short time scale which is associated with external probe- 
hadron interaction, and the other is a long time scale asso- 
ciated with the bound state nature of the hadron. In the QCD 
asymptotic region the short time scale is emphasized. However, 
till now, there is no prescription as how to incorporate the 
long time effects in QCli processes. 

There is another difficulty of infrared divergence and 

3 

mass singularity in QCl hard processes involving quarks and 
gluons. Quarks and gluons are supposed to be massless. 

Infrared divergence appears because of low energy massless 
gluons in the final state and the mass singularity appears 
in the quark propagator. Recently it has been showif that 
these singularities factor out and hence can be absorbed in 
the quark and gluon distribution function in hadrons. Many 
Calculations have been done involving hadrons in the asymptotic 
region of QCD and there seems to be fairly good agreement 
between theoretical calculations and the experimental data. 

We mention a few of them. 

(i) Deep inelastic lepton-hadron scattering ; The main 
conclusion of deep inelastic lepton-nucleon scattering experi- 
ments is that the nucleon structure functions i'2 show 

R o 

Bjorken scaling,^ i.e. ^2^^^ 

i'^(x,Q2) Fj^(x) 


2 

as Q o» , 
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2 2 

where Q = -q. - momentum transfer squared 


X = 


Q 


2P'q 


— r j P = 4--momentum of the nucleon target. 


But from the asymptotic freedom of QCD one knows that at any 
finite Q , however large it may be> some interactions among 
quarks are always present and hence one expects some scaling 

violation in deep inelastic lepton-nucleon scattering. Using 

6 

QCD the order of scaling violation has been estimated and it 

7 

agrees with the experimental data. 


e'^e" .Annihilation ■, The ratio R = 

cr( e e'-^Y P-”) 

O 

has been computed taking leading and next to leading QCB 
corrections and comes out to be 


R = 3 


f 


[1 


-f 


71 


+ (7.35 - 0.442f) ( 


a (Q' 
% 


2 

] 

(1.5) 


= 3.84 for 4 flavours, 

9 

which is in good agreement with the experimental value of 
4.17 + 0.09 (stat) +0.42 (syst). The factor 3 on the right 
hand side of (1.5) is due to the fact that each quark has 
three colours. 


(iii) HX ; Hadron-hadron hard scattering producing 

particles with large transverse momentum q^, ; In the framework 
of QCD these hard scattering phenomena arise due to quark-quark, 
quark-antiquark, quark-gluon, anti quark-gluon and gluon-gluon 
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hard elastic scatterings, the scattered q^uarks and gluons frag- 
menting into hadrons with large Gluck-Owen and Reya^^ have 

calculated the single pion inclusive cross-sections in p-p, 
p-p, TC-^ collisions in QCD hard scattering formulation using 
scale breaking Q -dependent quark and gluon distribution and 
fragmentation functions. The theoretical curve of differential 

t 

cross-section against q^ fits the observed ciurve. It is worth 
noting that the scale breaking effects play a dominant rol«. in 
these calculations. 

(iv) Jets : In any hard scattering process the predominant 
QCD radiative corrections are those due to collinear glue- 
strahlung and paircreation . These give rise to jets of 
partons. Thus two jets in e'^e’’ -*- y*-*- three jets in e’'’e”-^ 
Y* -►llg and so on are expected. Production of an additional 
jet costs an extra factor of a^(Q.‘^)i hence in e e"" annihilation 

o 

one expects 

P 

cr( 2 jets)*o(3 jets)io(4 jets) = l:a„;a^ 

o S 

For 3 jets in e'*'e~ annihilation, it is expected that there 
should be a large q^, cross-section. Recently evidence for 
three jets has been found in e'*’e~ annihilation,^^ which gives 
another credibility to QCD. ' 

(v) Massive dilepton production in hadronic collisions : 

p p X. This is another process in which the role 

12 

of QCD can be tested. In the Drell-Yan picture a quark 


from one hadron aimihilates an antiquark from the other 
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hadron to produce a massive photon which decays into 
lepton pair qq, y* First order QCD effects such as 

^q Kf + - lY* + - afe found to be of equal 

® -►p p ■►p p 

import ance.^^ 

In this thesis we have investigated some other aspects 

of dilepton production in the QCD framework. In Chapter 2. 

we give a survey of experimental status of dilepton production 

12 

and check it against the theoretical Drell-Yan predictions. 

In Chapter 3 we investigate the angular distribution of the 

dilepton in their rest-frame and find that it differs signi- 

2 

ficantly from the Drell-Yan prediction of 1 + cos 8 dependence. 
In Chapter 4 we study the asymmetry associated with dilepton 
production in longitudinally polarized hadronic collisions 
(pp and pp collisions). The asymmetry is quite large and 
varies with q^i, the transverse momentum of the dilepton, which 
Can be compared with future data and in turn may throw light 
on the performance of QCD in dilepton production. In Chapter 5 
we derive relations based on QCD expressions between dilepton 
cross-sect ions in nucleon-meson collisions, which in turn may 
be used to determine the distribution functions of nonstrange 
valence quarks and gluons in it and K mesons without requiring 
the knowledge of sea and gluon distributions in nucleons. 
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Chapter 2 


DILgPTON PRODUCTION IN HAPROIIC COLLI SIGNS 

In this chapter a short review of experimental situation 
of massive dilepton production in high energy hadronic colli- 
sions 

-► \x\~ + anything (2-1) 

will be given. We will also be comparing the experimental 

12 

results with the predictions of the Drell-Yan model in which 
an onshell quark from one hadron annihilates with an onshell 
antiquark from the other hadron to produce a virtual photon 
which eventually decays into a lepton pair, (Pig. 2.1) 

Assuming that the quarks do not have any intrinsic transverse 
momentum, the Prell-Yan differential cross-section is given by 

^2 — do ^ 1 . N ^ ^ cos^e) I e.^Jdx, f dx 

dm'^d cose 2m i 

„2 

[f^(x^) fT - (Xg) + ^{^2.^2 ~ s ^ (2.2) 

where 0 is the polar angle of one of the muons measured from 

the beam direction, in the t channel dimuon rest-frame, ^j_(x) 

is the probability of finding Che ith quark inside the hadron 

2 

with momentum fraction x and is the fractional charge, m 
the invariant mass square of the dimuon and /fs the c.m. 




9 
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energy of the colliding particles. Some of the predictions 
of this simple model ares 

(i) the differential cross-section is a fmction of 

m2 

single variable % = i.e. it scales as x. 

s 

(ii) the angular distribution of the muon is of the 

2 

form 1 + cos 9. 

(iii) the dimuon carry no transverse momentum. On 
introduction of intrinsic transverse momentum '\'600 Me7, part 
of which may have its origin in the fermi motion of the quarks 
inside the hadron, the maximum transverse momentum the dimuon 
can carry will be of the order 600 MeV in the Drell-Tan picture. 

Now one has to check these theoretical Drell-Yan predic- 
tions as well as the overall differential cross-section 
against the experimental data. Data on proton-nucleus colli- 
sion and pion and nucleus collisions at various energies over 
a range of dimuon mass are available. We enlist the salient 
points of the data as follows; 

Pro ton -Nucleus collision : Data^^ on the inclusive production 
of massive lepton pair in the mass range 2 Gev < m < 15 Gev 
by 200, 300 , 400 , 2000 Gev/c proton on nuclear targets are 
available at rapidities y = 0, 0.23 to 0.6. The following 
points emerge out of the data: 

(i) Scaling seems to be valid to the 20^ level. 
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(ii) The average transverse momentum <q.gi> of the dimuon is 
independent of mass hut rises with beam energy, Por 200, 400 
and 2000 Gev/c protons <qQi> is respectively 1.01, 1.2 and 
1.8 Gev. The data seems to be consistent with the fit 


<q^> = 0.6 + .0 22V s 

^ — 0.7 + .0018 s 


B 


d^q 


q.. 


T ^2^-6 


= C [1 + (^)^] 


(2.3) 


(iii) The angular distribution of file muon is being 

measured in the rest frame of the dimuon, the reference axis 

for 9 is either the beam axis (t channel Gottfried- Jackson 

15 

frame) in which we denote it by 6, or the Gollins-Soper 

axis which averages the beam and the reverse target momentum 

directions, in which we denote it by 6 . The Collins- 

c . s • 

Soper axis is suited for the case when one is incorporating 
the intrinsic transverse momentum of the partons also. 
Experimentally the two ^axes do not differ much. In the t 
channel frame, where we are ignoring the intrinsic transverse 
momentum, do /d cos 9 is of the form 1 + cos 9, whereas in the 
Collins-Soper frame do/d cos 9 is of the foim 

1 + ~ + (1 - I A^) cos^ ®C.S. 

P P 

% 2<K^ >/m , being the intrinsic transverse lODmentum of 
the parton. At beam energy 400 Gev/c and m ^#4-10 Gev 

<q2i> vl.2 Gev, kjp will be 'vl Gev. Hence •! fo .02 and 

2 

(dq/d cos 9„ ) *^1 + cos 9 „ • 
cs cs 
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According to Antreasyan, the angular distribution 
fit at y s = 6'2 Gev is of the form: 

1 + (1.6 + 0.7) cos^ 9 for 6 < m < 8 Gev 

c s 

1 + (0.5 +0.6)- CO 0 for 8 < m < 11 Gev. 

M. C S 

Though the authors claim that above the threshold the 
angular distribution seems to support the Drell-Yan picture, 
yet the picture is not very clear because of large uncertainty. 

(iv) Shape and norm.alization of the differential , cross-section 

: It all depends on the structure of the sea distribution 
dm 

in nucleons. If one assumes a symmetric sea u(x) = d( x) = 

s(x) = s(x) = c(l - x)^/x where p = 7 or 9 ard c is fixed by 

the momentum carried by the sea, then the shape of agrees 

dm"^ 

with the experimental curve, over all normalization being 
arbitrary. Overall normalization will be all right if the 
sea distributions extracted from pK-^ p’^p'x data and vH-^-px 
vN -►laZ data are identical at the same values of x and Q (=m ). 
The neutrino data for 0.05 < x < 0.25 and 5 £ £ 20 Gev^ 

overlaps that of the pN p’^p'^x data only for x near 0.2 and 

2 2 2 16 
Q'^(=m ) = 20 Gev . E. 1. Berger determined the magnitude 

of sea from the two data. In the overlap region, the magnitude 

of the sea determined from v If data is at most one half of 

that determined from the pN-^ data. 
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(v) Tal’get dependence: It is assumed that in the asymptotic 

region partons interact incoherently, hence it is expected that 

the dimuon cross-section should vary linearly with A, the atomic 

number of the target. The data^^ indicates that the target 

2/3 

dependence changes from A ' dependence for m£ 4 Gev to an A 
dependence . 

14 

Pion-nucleon collisions : Data on the inclusive production 
of massive dimuons in the mass range 2 Gev < m < 12 Gev by 
150, 175, 200, 225, 280 Gev/c pion beams on nuclear targets are 
available. The salient features of the data are given below: 

(i) With the exception of 225 Gev/c data of Anderson et al^^ 
scaling is well satisfied within the normalisation error of 
+ 15 ^. 


(ii) For dimuon mass m< 4 Gev, the average transverse momentum 
<q.Q,> seems to rise with m. Above 4 Gev, becomes indepen- 

dent of dimuon mass but it increases with increasing beam 
energy. For instance at 150 Gev/c beam energy it is 1.04 Gev/c 
whereas for 280 Gev/c it is 1.2 Gev/c. The data seems to 
give the following fit 


<Ggi>=a + bVs a».7 Gev/c 

b = .02 


1_ dg 



F 


D = 0.49 + *05 nb /Gev/c 
S = 1.7 + .05 Gev/c 
F = -3.2 + 1.3 
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= A [1 - IXj - B|] = 

A = .43 + .03 nb 

B = .14 + .02 

C = 2.1 + .03 (2.4) 

(iii) In the pion-nucleus collision, the valence-valence anni- 

-f 4- — 

htlation is the dominant process. Hence in n H.*- p p. X 

annihilation ‘gives the important contribution 
and in ti-H ♦p‘^p’'X it is the '^1(2^2) annihilation gives 

the major contribution, and according to the Drell-Tan formula, 

• 4. 4. 1 4.1 4. r. u'^u'’X) . ®d 1 

it IS expected that H = — ^ — 2 ► — ~ 4* 

o('n:”H-*' p p~X) large ^ ^ 

14 

, 1‘he data of Anderson et al confirms this tendency: for 
low dimuon mass m, E 'vl and near m '\<8 Gev it is 

(iv) The angular distribution j-”- channel rest frame’ 

2 

of the dimuon is found to be of the form 1 + a cos 9. The 
value of a as fitted by different groups is not very consistent.; 
According to the- fit of Badier et a = 0.80 + .17} Barate : 

et al.^^ give a = 0.52 + .46 but according to inderson et al^"^ 
a = 1 at low X but a = -1 near x ^ul. The last quoted result 
gives a clear evidence that the dimuon production in x 'v 1 1 

region is not predominantly through on-shell quark annihilation i 

(v) To check the absolute normalisation one extracts the i 

valence quark ( antiquark) distribution in pion from -jxH .»■ p’*p”X 
data using the quark distribution in nucleon determined from ; 
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deep inelastic scattering and then checks with the established 
fact that of the pion’s momentum is carried by valence quark; 

Data of various groups are not consistent, hewmann et al.^'^give t; 
fit f^(x) = a(1-x)P/x, p = l.ai + 0.05, A = 0.52 + .05 for 
<in^ =! 25 Gev^ and find / x dx « 0.2, showing that 40^ 

of the pion's momentum is carried by valence quarks and anti- 
quarks which seems to be reasonable. I'he later data of Barate 
14^ 

et al . gives / f (x)dx = 2.8 which is 2.5 times higher than 

0 ^ 

expected - which shows a clear disagreement with Drell-Yan 
model. If higher order QCD effects can increase the cross- 
section .by 100% the result can be consistent. Badier et al.^^ 
also report that the experimental dimuon differential cross- 
section is almost twice of the Drell-Yan prediction. 

17 

(vi) Target dependence: The Chi cage -Princeton collaborators 
find an exponent a = 1.12 + 0.05 for m > 4 Gev, whereas recent 
data of Badier et al^^ supports the idea of incoherent parton 
interactions with a = 1.03 + .03. 

Thus we get the following important information from the data 

■"1“ M. — 1“ 

on pN-»- p p X and Ttiff p p x: i 

(1) that the scaling seems to be valid, i 

(2) above m = 4 Gev, the average transverse momentum of the | 

dimuon <Q.2i> independent of mass but increases with fs - a 
fact unexplainable by simple Drell-Yan model. | 
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(3) the angular distribution — ~ is of the form 1 + acos^e 

CL 0000- 

where at low x, a seems to be equal to 1 but in the region 
X ''' 1, a = -1 indicating a clear departure from the simple 
Drell-Yan model. 


(4) Ihe sea distribution extracted from pN-^ p \i~ X data is 

1 

twice of that extracted from v JT -+■ jiZ data, / f (x)dx extracted 

+ - ° 
from tcN-*- p p X is 2.5 times greater than the expected value of 1. 


It is clear that the Drell-Yan model gives only a quali- 
tative but not a quantitative picture of the dilepton production 
in hadronic collisions. Leaving aside the absolute normalisa- 
tion, attempts have been made to explain the high transverse momentu 
momentum of the dilepton by introducing intrinsic transverse 
for the quarks. Of course the quarks do have a transverse 
momentum distribution because of fermi motion and quark-gluon 
interactions inside the hadron. If one assumes a Gaussian type 
of distribution for quark intrinsic transverse momentum as I 


fCS-p ) = e 


~ 2 
o/ 




( 2 . 5 ) 


then from the Drell-Yan formula, one obtains 


d^a 


2 2 

dm^d^qj. 




■dm‘ 


<q, 


T 


> = 


a. 


= < 


% 


>, <q2.> 


% ° 

2 o 


( 2 . 6 ) 
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which is decreasing too fast to reproduce the experimental 
curve (fig. 2.2). Moreover, one reauires an unlikely large 
average intrinsic transverse momentum < Kj> 'v 1 Gev/c to explain 
the observed <qQi> "v 1.4 Gev/c, and this does not explain the 
increase of <q.,j, > with incident energy. 

Now it is quite natural to look for the higher order 
QCD effects in the production of massive lepton pair. To the 
first order of QCI) perturbation, either a quark from one hadron 
Can annihilate with an antiquark from the other hadron to 
produce a gluon and a virtual photon decaying into a lepton 
pair (qq-t-qy fig. 2.3(a)) or a quark (or antiquark) can 
interact with a gluon from the other hadron to produce a quark 
( anti quark) and a virtual photon decaying into a lepton pair 

fig. 2, 3(b)). For processes where valence-sea contri- 
bution is important such as pp collision the second order 
quark-quark scattering (qq-^qqy*) producing a virtual photon 
may be comparable to the first order effect because of the 
dominance of valence quarks. These QCD processes can give rise 
to the transverse momentum of the dimuon because the final 
state recoiling gluon (in the case of qq QCD annihilation) or 
quark (in the case of qg compton process) can balance the 
transverse momentum of the dimuon and thus can explain, after 
taking the intrinsic transverse momentum of the part on into 
account, the observed transverse momentum of the dimuon. 
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Explicit calculations on QCD contribution to inclusive dimuon 

13 

production in hadronic collisions by Piitzschand Mnkowski, 

13 

AltafSlli, PaTisi and Petronzio and others show that the 
quark gluon QCD compton process in pp collision and quark- 
antiquark QCD annihilation in the pp or Tip collision is as 
important as the Drell-Yan process and the average transverse 
momentum <qji> of the dimuon due to these QCD processes 
increases with Vs at a fixed x. Ihe transverse momentum 
distribution of dimuon due to first order QCD effect is given by 


a(s,m^,q|) = — ' a = T e.- //dx.dx e ((g _m^) ^ - 4sqm^) 

dm^dqn, i ^ ^ ^ 


d^0 


2 n 2 2 . 






{ fi (xj f^ (xj + (i^i)} + -Mil- { 


dm^dq/ 


2 2 
dm^dqj,^ 


(fg^(x^) + % (2:2) + ^^2^ 


(2.7) 


where 

. 2 .-.„ 2 


l6a^a 


'‘2 4 

( ! -o- - 2) 


dm^dq^^ 27 % V( s-m^) 2-4sq^^ 


m 


2 d"o 


2;^g 


2 

a a 


.2,„ 2 


dm dqr^, V( s-m^) ^-4sq,^' 

(1 - 


rs + 3m^ , (s-m^) 

L — + 9 


and s = x^x^s . 


( 2 . 8 ) 
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'The QCD contribution is not enough to produce the observed 

q.^ spectrum of the dimuon - the contribution due to the intrinsic 

9 

transverse mom-entum of the partons is aJLso needed. The intrinsic 

P P P P 

contribution is added to a(s,m,q,j ) after convoluting a(s,m ) 

with the quark premordial motion which is assumed to be 

f(K^^) = — — g e ^ xfhere 

27ra^ 

Thus one obtains 


total 


(s,m2,q^2^ ^ f(q„2) + / d^i^ a(s , m^, ( ) 2) f ( 

dm 


(2.9) 

Taking = .68 Gev <Kj,> = = 600 Mev, the experimental 

19 4> — 

data on pp -»-p p X at ys = 27.4 Gev is exactly reproduced. 

As the QCD quarkgluon contribution in pH collisions and QCB 
quark- anti quark annihilation contribution in icN collisions are 


of the order of Drell-Yan contribution, the problem of normaliza- 
tion of Sea distribution in pH p'*"p~X data and valence quark 
(antiquark) distribution in pions in ttK’^ p''’p~X data may be 
taken care of. 


It looks as if QCD is playing a very important role in 
the production of dilepton and hence it is desirable to study 
some other properties of the dimuon production which may reveal 
something more about the QCD effect in dilepton production. We 
study the angular distribution of dimuon coming via first order 
QCD effects. At high transverse momentum the intrinsic momentum 
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of the partons can be forgotten and hence the angular distri- 
bution will be only due to first order QCD processes. Hence 
it is desirable to study the dimuon angular distribution at 
large in order to understand the QCD role in dimuon produc- 
tion. ¥e also study the asymmetry associated in the production 
of dimuon with polarized beam and target. I'he asymmetry can 
arise due to subprocesses involving particles with spin, the 
nonzero asymmetry will rule out the models of subprocesses 
involving mesons such as quark-meson scattering etc. Moreover, 
this study may lead to the knowledge of spin dependent distri- 
bution functions of partons inside polarized hadrons. 



Chapter 3 


U^CULAH DICTEIBUTION op dileptons PSODUCED in 
hadronic collisions 


3.1 Introduction 

12 

The mechanism of Dr ell and Tan for lepton-pair 
production in hadronic collisions, involving on-shell quark- 
antiquark annihilation (Pig. 2.1) into a massive virtual 
photon, is not adequate to explain all the details of the 
massive dilepton data, particularly the high transverse 
momentum of the dilepton pair. 

The first order QCD effects , viz., (i) qq -►gy'^ 

(Pig. 2.3(a)) and (ii) qg-»- qy* (Pig* 2.3(b)) in the produc- 
tion of massive dilepton have been investigated by Pritzsch 

13 

and Minkowski, Petronzio and others, and these first order 
effects in addition to the contribution of primordial trans- 
verse momentum of the partons give the observed transverse 
momentum of the dilepton and its variation with incident 
energy i.e. <q2i> = a + bfs. It is natural then to investi- 
gate some other aspects of the dilepton production which may 
give additional information about the mechanism of dilepton 
production, especially regarding the QCD contribution to 
the dilepton production. One such aspect would be the 
study of angular distribution of one of the dileptons in 
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the dilepton rest frame. In contrast with the simple Drell- 

p 

Yan angular distribution of the form 1 + cos e, one expects 
a nontrivial 6, ({> dependence in the angular distribution if 
first order QCD effects gy" and qg.^ qy* are taken into 

account and then consequently one can find out the asymmetry 
associated with the niomber of events on two opposite sides of 
a chosen plane, which can be checked with the experimental 
data. 

In this chapter we study the angular distribution of 

4 " — — 4 - 

muons in Hp p p X, H = p,p,m ,n by taking into account the 

first order QCD effects qq-*- gy*' and qg-^ qy'^ . ¥e calculate 

the coefficients A, B, C occurring in the angular distribution 

of muons in the p'^p” centre -of -mass frame, which is given by 

1 + A cos 9 + B sin 2 0 cos4» + C sin 0 cos 24, where 9 , 4 are 

the polar and azimuthal angles of a muon with respect to some 

specified axis. A, B, C can be functions of various kinematic 

2 

variables like s, m , q^, where q^ is the invariant generaliza- 
tion of the c.m. transverse momentum, aiid the prediction of 
their dependence on the kinematic variables, as well as their 
actual magnitudes can be tested experimentally. 

We obtain the coefficients A, B and C occurring in 

j 4 — o o 

— 2 1 as functions of s, m and q^ . We restrict 

dm^dqgi'^d cos0d4 

ourselves to large q^, (q^ > 1 Gev) for two reasons. One 
reason is that we neglect parton transverse momenta *^.6 Gev, 
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which via the Drell-Yan process would contribute at low 
For large qip, where the Drell-Yan process does not contribute 
significantly, the values of A, B, C become independent of 

p 

the QCD running coupling const anx a = g /471:, and so 

s s 

arbitrariness due to the choice of a is avoided. I'his is 

s 

2 

fortunate, since for time like q (which is the present case), 

the value for has to be obtained as an extrapolation of the 

value measured in inelastic scattering for spacelike q^. Ihe 

second reason is that by restricting to large q^, our results 

po 

are independent of mass-eingularity regularization procedure, 
whose effect would be important for low values of q^. Our 
calculations of angular distribution are done in the Gottfried- 
Jackson frame, i.e., in the c.m. frame, we measirre the 

angular distribution with the beam direction as the z axis and 
with the X axis in the plane of the two hadron momenta. 

21 

Kajantie, lindfors and Raitio have also studied the 

QCD angular correlations for dimuons with rapidity y = 0 . Our 
22 

study is for y-integrated distributions. Also, we have 
calculated the coefficient, we call B, in addition to A which 
they calculate. 

3 .II Derivation of .jAngular Distribution 

To study the angular distribution of a muon in the 

process 


%^^1^ H 2 (P 2 ) + anything (3.1) 
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(the quantities in the parantheses denote the respective four 
momenta of the particles), we must obtain the differential 
cross-section d o/dm dq^j-, d cos9d(j) for this process, 6 , <j> being 
the polar and azimuthal angles of one muon with respect to 
some specified axes. ¥e shall be working in the c.m. 

frame, and we choose the direction of as the z axis and 
the direction of P^ x P^ as the y axis. This is the Gottfried- 
Jackson or the t-channel helicity frame. Assuming the muonpair 
to arise from a virtual photon of mass m, we can conveniently 
express this differential cross-section in terms of the density 
matrix elements. Por the process a+b^c+y* ^ > the 

■*-\i \i~ 

matrix element can be expressed as 


M = a) ( a+b c+y"''' ( helicity=i) • ^ (y'''‘ ( helicity=i) 


which gives 

I |m| 2 = . 1 |i = p 

all spins q^ 

i,J — 1,0, —1 . 


(3.2) 

( 3 . 3 ) 


The matrix p in the restframe of the virtual photon of 4 
momentum q is given by 


2 2 
P = e^q"' 


^ 1 + cos ^9 


sin29 

e 


T2 


sin^e 


sin 2 6 „i|) 
V2 

2 

2 sin 6 


sin^0 e^^^ 


sinpQ 

— -" T T ' T — e ^ 


12 


sin29 2^ 

-^2 — ® 1 + cos 9 


( 3 . 4 ) 
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So the differential cross-section ( — S -- — ~) /{ — 

dm^dqr^.'^dn^ dm^dq^.'^ 

can he expressed in terms of the density matrix elements as 


’"''dm^dq^S > = L ["ll +<>00 

4 ^2f 


sin^ S 


+ f2 Re p^Q sinSQcos <j) + sin^9cos2<J> ] (3.5) 


where the density matrix is normalised so that 2p^^ + Pqq = 1. 
¥e shall henceforth work with unnormalised P. •, in terms of 

1 j . 

which, the distribution is 


1(9 , (})) = 1 + A cos^9 + B sin29cos(j) + C sin^9cos2(|> 
with 

^ ^ Pfl ~ PQO 
Pll PQO 

B ^ . I L ^10 . - 

^11 ^ ^00 

^1 -1 

Q ^ _ 

Pll ^ ^00 


(3.6) 


(3.7) 


The density matrix elements are sums of contributions 
from all incoherent processes, and they are integrals over 
panton distributions of the corresponding density matrix 
elements (denoted by p. .) for the subprocesses involving 

X tj 

partons . 

p for a particular subprocess involving two partons 
in the initial state can be written as 
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Pi j 


I 


(3.8) 


q“ flux “ ^ 

where is the amplitude for the production of photon of mass 
m, helicity i, aud fixed transverse momentum in the parton 
c.m. frame? fhe summation refers to averaging over all initial 
spins, summation over spins and integration over momenta of 
all final particles but the photons? and summation over initial 
and final colours and averaging over initial colour. Thus, 


fp 


11 


lii d^a/dm^dq^,^ 


(3.9) 


¥e shall be evaluating p. . for the first order QCD sub- 
processes (i) qq-*- gY"'" and (ii) qg-^qy'^ iii addition to the 
Drell-Tan qq annihilation process (qq-^Y'"')* 


The Drell-Yan process % The amplitude for the subprocess 

q(p) + q(k) Y"(l) (3.10) 

given in fig. 2.1 is 

® ®q ^i (3.11) 

where e^(q) is the photon polarisation vector for helicity i, 
and e^ is the quark charge in the units of the positron 

q 

charge e. Taking 

e^l = p (0, +1, -i, 0) 


( 0 , 0 , 0 , 1 ) 


( 3 . 12 ) 
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as the representation, for the polarisation -vector of a massi-ve 
photon a.t rest, we obtain 


.vDl 

hi 


2 2^- 

e a s ^ „ 

— 6(3-ia^) 6(qrf.^) = 


P -1 , -1 


(3.1: 


/s DY 

remaining p . . are zero. 

J_ J 

■where the factor i comes because summing and averaging over the 

^ 1 

initial colours of the quark and antiquark i.e. T 1/3x3 = v , 

2 2 2 ^ 

Vs/2 = part on c.m. energy and m"^ is the in-variant 

mass squared of the photon. 


It is clear that for laxge qrj,? with which we shall be 

DY 

dealing, the contribution of p will be zero. 


The subprocess qq.>- g-f' 

The amplitude for the subprocess 


q(p) q(k) ->• g(p’) + Y'"(l) 
according to the diagram (fig. 2. 3a) is given by 


(3.14) 



Ss® 



e*"^(p') 


eV (q) v(k) 


+ X 
X ”X 


^ p-p' ^ 


+ Y 
V 


^ 

P'-P 


Y,,]u(p) 

1-1 


(3.15) 


where g is the strong coupling constant, £®^(p') is the polariza 
tion vector for the coloured gluon a (a =1,2, ...8), X®" are the 
usual Gell-Mann matrices for SU( 3) normalisation 


Tr[x^X^] = 26 ^^ 


(3.16) 
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and X is the SU( 3) colour spinor. Using (3.12) for eV(q,) and 
summing over all • gluon final states we can express in terms 
of the folio v/ing invariants: 


s = (P3_ + - 2P£P2 

2 2 

q = m = invariant mass square of the virtual photon 
s = (p+k) = (x3_P2+X2p2) “ px^x^s 


t = (p-p’)^ 


u 


(q - p) 


^ ^ A. p 

3 4-1; + u = m 


M 


q^ = V”^ = Transverse mom. of the virtual photon 


V S 


= 


21 

Ys 


X = m^/s 

x^jX^ are the fractions of the hadron momenta 
carried by the partons. 


(3-17) 


We get the following expressions for p'?’^ 

1 3 


'qq 

^00 


;^q 

^1 


= 2 


2 2 / 4 ^>2 2\ 
, a a„ (m +s - 2s qm ) 


Pl-1 


2m^ 


-^qq _ q " "s x ■ r 

- ^ 2^2 ^ 1^11 ^iii 2 2\ 

xSte ms V ( s*"2! ; -4sq2-i ( m + q^ ) 


(3.18) 


m 


2 2 / 4 -"■2 2\ 

.qq _ eg a a3 (m +s -2s q^ ) ^ 3^ 

IStud^^ V’(s-m2)^-4sq„^ (m^+q|p^)^ 


] (3.19) 


/^qq 

PlO 


(m^+s^-2s q^^) ^ m(q;p^ - m^) 


PlO^" = — ^ 7-2- - T - 2 

V’2qrn(m +qm ) 


] ( 3 . 20 ) 


1871 m s'^f ( s-m^) -4sq 


•T 
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The above expressions have been obtained on the assimiption that 
the quark momentum p is along the z axis, that is, the quark 
comes from the hadron whose beam direction is along the z axis. 
However, it turns out that the contribution when the quark 
comes from the other hadron and the antiquark q is along the 
z direction, is the same as i.e.. 


Here a 
All the 

— y 

5x3 } y, 

Q, y D 


s 


pM = pM 

~ Sg q^^ ^ 


expressions for p 




C - H 


(s-m^)/2s . 

include the colour factor 
2x8 4 

_ 


(3.21) 


The subprocess qg'*‘qY' ' 

The amplitude for the subprocess 


q(p) +g(k)->q(p') +-f(q), 


(3.22) 


obtained from the diagrams in Jig. 2.3(b), 

= e e g £t'^(q) £^’^(k) u(p') [y Y, 

^ q s 3 . p > 

, Xa 

u(p) [X ~X] 


+ Y 


V 


_1 



(3.23) 


Notations are the same as defined in the previous section. 
Using (3.12) for e^(q) and integrating over the quark final 
states in (3.8), we get for the following expressions: 
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;q.g ^ 2^"^^ 


00 


1,-1 


^22 


967rs^m^f ( s-m^) ^-4sq2i^ 


.„2/ , 2\ 2 

4bi ( s+m ) 

(2 ^ „ 2 n 2 

(m + ; 


4m ^( 3s+m^) 

2 2 

m +qj, 

(3-24) 


:ig _ 


ai 


= p 


1,-1 


= 22 


96718 ^m^y(s-m^) ^-4sqj 


r 8 - „2\ 2 . „4i 

[- _t ( s-m ; + m / 


2 S q„ 


+ ( 3m2+s) ] - I 


(3-25) 





[_2(s+m2) 

2qrj- 


o s + m^) ^ ( 8^ + 8sm^ + 3m^) 

(m^+qrj,^)^ (m^ + q^,^) 


(3.26) 


The above expreseions have been obtained on the assumption 
that the quark momentum p is along the z direction, that is, 
the quark comes from the hadron whose direction is defined as 
the z axis. However, we have also to consider the contribu- 
tion when the quark comes from the other hadron i.e. the gluon 
momentum "S; is along the z axis. ¥e obtain the following 
expressions for the corresponding density matrix 


^gq 
b 00 



r 8m^( s+3m^) 

^ 2^„ 2 
m +q^ 


8m^( 3+m^) ^ -1 
s(m^+q^^)^ 


(3.27) 



< CL 
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^ gq _ gq 

Pii " p -1,-1 - 


a a. 


96txs^iii^^( s-m^) ^-4-s ® 


^ 2 

S— m f / ^ 

[- — _{ (s_ni ) 


^ „ 2^2 


+ } 


+ (^Ul-rS>J -^-p 


gq 

oo 


("3 


gq _ gq y:;- _ 
10 - ■^p-1,0^ ~ 



^ IQ \ ^ r 8^ 1 S L 

2q^^ s(m^+qj^)^ 



4- 2(3s + 5ni^)] 


(3.29) 


All expressions for and include the colour factor 

■ 5 ^ 1 X ^/2 ^^/ 2 ) = ^ which arises because we sum over 

s , b 

the colours of initial and final quarks and gluons and 
average over the initial colours. 


To obtain the complete density matrix elements for 

the process H^E 2 -»- y'^X, the have to be convoluted with the 

appropriate structure functions and then added together. 

If we denote by f|(x) and f|(x), gj^Cx) the structure func- 
tions of a quark q, an antiquark q and a gluon respectively 
inside the hadron H, then the convo^-uted density matrix 
elements are given by; 



I dx. 


1 

J dX 2 9 


(( 


Tn2\2 

-m ) - 


■4s qj'^) 


X J 


"g^(Xi) 


f| (Xp) 


(3.30) 
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q.q 

Pi3 


^13 


gq 

Pid 


1 1 


= / J ((s-m^)^-4s q^.^) f| (x^) fg {x^) 


0 o 


1- 3 


(3.31) 


fdx^ Jdx^e ((s-m^)^-4s 


^1 J v.>u^p 

0 o 


(3.32) 


1 1 

Jdx^ Jdx^e ((s-m^)^~4% q^^) p^^ % (^q) ^H ^^2^ 


o o 


(3.33) 


,<lg ^ 
id 


go. 

Pid 


11^ 

/ dx^ / dx^d ((s-m^)^-4s 1^^) P J ^ 


(3.34) 


1 1 
Jdx 

0 0 


Jdx^ Jdx^ 0 ( ( s-m^)^-4s 11^)^^^ 


"1 


(3.35) 


The last two expressions (3-34) and (3.35) are the contribu- 
tions of anti quark-gluon scattering. The Crell-Yan contribu- 

2 

tion is zero because q^-, >0. Hence the density matrix 

elements p^ ^ for the process Y*X at q^ > 0 is 


p. . = I (P ^<1 H-P^'l +P^^' +P +p'^g + pS*^) 

aai flavours 

(3.36) 


. . are functions of s, and which may he expressed in 
1-3 « !• 

terms of s, t = 


and r„ = 




T " fs' 


Once this is done the 


s dependence is factored out and is the same for all p. - 

X y 

Knowing p. A, B, C can be computed using equation (3.7). 
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A, B and 0 are functiom of the sealed variables t = m^/s 

ard r^ = q^i/Vs, because the s dependence and a^{n^) 

cancel out in the numerator and denominator of eq. (3-7). 

Of course, we are treating the structure functions to be 
2 

q -independent . 


3. Ill Numerical Calculations and Results 


¥e have- evaluated A, B and C from p. given by 

1 U 

equation (3.36) using (3.7) for various values of t and r^ 
for the processes pH -*■ H = pjp, tc”. In our 

calculations, we have used the expression 


agCm^) = [1 + 


“o 


( 3 .? 7 ) 


given by QCD and the renormalization group equation for three 

colours and four flavours. We take ag(mQ)=0.25 for mo=2 Gev. As 

•we restrict ourselves to nonzero transverse momentum of the 

dimuon, the Drell-Yan process does not contribute to P. . 

1 '«j 

2 

and consequently in the calculation of A, B, G, o:_(iQ ) 

s 

cancels out and our results become independent of the 

o P P 

choice of a (m ). We have neglected the m (=q ) dependence 


of the structure functions. 


I he process pp'» ; Por evaluating for this process 

we need the quark and gluon distributions in protons. Por 

the quark distributions, we use the parametrization of 

23 

Peierls, frueman and Wang : 



56 



11 

= 1.79 

(l-x)^(l+2.3x)/fx + s(x) 

F(.) 

= dp(x) 

= 1.07 

(l-x)^'l/Vx + s(x) 


= 

= fp^-) 

= s (x) = 0.15 (l-x)'^/x (3.38) 

It 

For the 

gluon distribution, 

we use 


gp(x) = (l-x)’“/x (3.39) 


with m = 5, and 8. The factor — takes care of the fact 
that gluons carry 50^ of the proton momentum. 

The values of A and B as functions of r^ are pre- 
sented for various values of x in Tigs. 3.1(a), 3.2(a) for 

g (x) = 3(l-x)^/x and in Tigs. 3.1(h) and 3.2(b) for g (x) 

P P 

= 4-5(l-x)®/x. (C is given by C = ^(l-A).) A, B do not 

2 

separately depend on s, m and scaled variables 

T and r™. For evaluating p. we have used ][s = 27.4 G-ev. 

The shape and size of A and B do not change much for the 

5 8 

two different gluon distribution 3(l-x) /x and 4-5 (l-x) /x, 
although P. ^ change appreciably. A and B show interesting 
variations with r^^,. For low values of r^ <0.05, A is 
almost 1 and then decreases with r^ to a minimum of about 
0.1 and then increases. B starts with a value near zero 
for r^ = 0, decreases with increasing r^, to a minimum of 
about -0.12 and then increases to a value of about +0.2, 
after which, the curve flattens, and B decreases slowly. 




ri6, 3 . HQ) 
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These variations of A and B with r^, clearly indicate that 

(i) at low Pg, < .05, the angular distribution is almost 
01 the form 1 + cos 0 and hence low transverse momentum 
data on dimuon production is not fit for the investigation 
of the QCD effects. 

(ii) for r^ > .05, there is quite strong dependence on 
both 0 and (J) . This region is suitable for the investiga- 
tion of the QCD role in the dimuon production. 


A is measured in the <{> -integrated distribution, which 

p 

is of the form 1 + A cos 9. An interesting way of observing 
nonzero B is to consider the difference in the fractional 
number of events occurring in the two sides of the yz plane, 
with 6 restricted to values O<0 <9^0_< •n:/2. 

For 0 ^ = 0, the corresponding asymmetry is 


6(0^) 


©1 5it/2 

sin0d0{/ l(0,ii>)d<}>- J 

m/2 -m/2 

01 . 5m72 ; ; fj2 


J 


m^2 




slrB d0 { 


m/^ 


l{0,(J.)d(J. + J I(8,<j>)d4>} 
-m/2 


^oot 

m 


.2 
sin 0 


1 


^ 3+a(1+cos ©2^+ cos^e^) 

^ is, thus, proportional to B. For 0^ = m/2 


(3.40) 


4.B 

m( A+3) 


6( m/2) 


(3.41) 



42 


Por = 71/3 

6(tx/ 3) = _ 1 2+f a) (3.42) 

At T =0.1 and = 0.3, for example, A = 0.18 and B = 0.12, 
6('ji;/2) - -.06, 6(n;/3) - -.07. Thus the asymmetry can be 
quite large. 

18 

As pointed by Georgi it is a bit dangerous to 

ignore the higher order QCD subprocess qq-^-qqy* + _ 

+ - \i-~ 

pp -*• p \i X, but fortunately, an explicit calculation by 

24 

Eripfganz and Contogouris shows that the contribution 
due to this subprocess is negligible over most of the 
kinematic region. 


+ — 


The process pp.»- p p X : To evaluate P . • for this process, 
we need the parton distributions in antiprotons in addition 
to the parton distributions in protons given by (3.38) and 
(3.39). The parton distributions in antiprotons are related 
to the parton distributions in protons by charge conjugation. 


f-(x) 

p 

= fp(-) 

= u ( x) 
p 



fl(x) 

= 

11 

pb 



f-(x) 

p 

= f“(3t) 

p 

= 4 m = 

f|(x) 

= f|(x) = f®(x) 


fp(x) = f®(x) = s^(x) 


(3.43) 
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The gluon distribution in p is assumed to be the same as 
that in p. 

gp(x) = gp(x) = ^(l-x)“^/x (3.44) 

We compute p. ^ using the functions u (x), d (x), s (x) 
IJ P P P 

given by equation (3.38) and then calculate A, B and C for 
various values of t and r^. The plot of A and B against r^ 
is shown in figs. 3.3(a) and fig. 3.4(a) for g (x) =3.0 

P 

(l-x)^/x. The shape and size of A and B do not change if 

O , 

we change the gluon distribution to g (x) = 4.5 (l-x) A- 

P 

This is understandable because for pp-+ the quark- 

gluon comp ton subprocess is negligible in comparison with the 
quark-ant iqu ark QGD-annihilation subprocess. Following 
comments are in order. 

(i) The shape of A, B versus r^ curves for pp-^p'^}j.“X are 

similar to those in pp process. 

(ii) A starts at -x/l for r^^ 0.02 and goes to a minimum of 

0.2 and then rises with r^. 

(iii) B s 0 for r^ ss 0 and then decreases to a minimum of 
'u- 0.2 and then increases with r^ and goes up to . 22 . 

(iv) The angular distribution for r^ < .02 is of the form 
1 + cos^ 0 , same as the Drell-Yan prediction. Hence this 
region is not suitable for the study of QCD effects in dimuon 
production. 
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(v) Por > 0.05, both 8 aiid <{) dependences are appreciable 
and quite different from 1 + cos^q behaviour and thus large 
r^ region is sui'ba’ble for the study of the role of QCD processes 
playing in the dimuon production. The asymmetry 6(9^) defined 
in (3.40) is appreciable for this process too. Por instance 
at T = .1 and r^, = 0.2, A = -22 and B = -.15 and hence 
6('n;/2) = .06 and 5i%/3) - .08. It is to be noted that the 
asymmetry 6(6^) can change sign if r^ is changed because B 
varies from -ve to +ve value as r^. increases. 

The process -n;'''p -» : Here, in addition to the part on 

distribution in proton, the Imowledge of parton distributions 
in pions is also required. ITo parametrization based on 
deep inelastic scattering on pions, due to the lack of data, 
is available. We use the following parametrization due to 
Parrar^^, derived on the basis of SU(2) symmetry, charge 
conjugation symmetry, quark counting rule and the fact that 
40^ of pions momentum is carried by valence quark and the 
rest by the gluon {50fo) and the sea 

f^+(x) = = f^I(X) = f^^(X) = u^(X) 

% % % % 

= 0-4 (l-x)/x for X 0*5 

= (0.1 + 0.72Vx)/x 





4a 




FIG. 3.6 (<J) 
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S^(x) = f^®(x) = = f'^f(x) = f'^^Cx) = f®J(x) 

% % % n 

= f®f(x) =0.1 (l-x)^/x 

71 

g ^(x) = g (x) = 2(l-x)^/x 

X 71 

u^(x) = uj[(x) + s^(x) (3.45) 

The superscripts v and s refer to -valence and sea quarks. 

Using the parton distributions in pion given by (5.45) 
and in proton given by (3.38), we calculate A, B for Tx^p+-|j,'^|i”X 
for various values of t and r^. Bor various values of t, a 
and B are plotted against r^. Pigs. 3.5(a), 3.6(a) show A and 

jr 

B against r^ for gp(x) = 3.0 (l~x), /x respectively. 

No appreciable change is noticed in going from 

{r Q 

gp = 3,0 (l-x)^/x to gp(x) = 4-5 (l-x)°/x. Following points 

are worth mentioning regarding the variations of A and B 
with r^,. 

(i) The shapes of the curves are similar to the curves for 

pp-*- 

(ii) Starting from'ul at r^,''* 0.05, A decreases to a minimum 
of ''-0.15 and then increases again with r^^. 

(iii) B is almost zero for r^ < .02 and then decreases with 
r^i to a minimum of -0.15 and increases and flattens 
beyond r^ ^ .3. 
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(iv) The angular- distribution for low r^ (<0.02) is of the 

2 

form 1 + cos 0, same as the Drell-Yan process. 

■'(v) Beyond r^, = 0.1, effects of QCD show up in the dimuon 
angular distribution, the angular distribution being 
quite different from 1 + cos 0 behaviour. The asymmetry 
6(0^) for T = .1 and r^ = .3 is 6(Tt/2)« -.03- The 
asymmetry will be large for low value of i:. 

The process 'n;~p-*- ; ¥e calculate A, B and C for this 

process for various values of x and r^ using the parton dis- 
tributions given by (3.38) and (3.45). In Pigs. 3.7(a), 3.8(a) 
for g^(3t) = 3(l-x)^/x, we show the variation of A and B with 

r^ keeping x fixed. There is no significant change if we take 
g (x) = 4.5(1 -x)®/x. 

Jr 

It is found that 

(i) the shape of the curves are similar to those in pp^^i'^p.'X 

(ii) A starts near 1 for .05 and goes to a minimum of 

0.2 and then increases with r^. 

(iii) B is''/0 near r^^. 'v .02 and takesa minimum value of ''•-0.2 
and then increases with r^ and goes up to +0.22. 

(iv) The angular distribution for low r^ (<0.02) is of the 

2 

form 1 + cos 6. 

(v) Por large values of r^ ( > -05), the angular distribution 

■ 2 

is different from the Drell-Yan prediction 1 + cos 0, 
showing a significant role being played by QGD effects 
in dimuon production. 
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(vi) The asymmetry 5(71/2) at r,^ = .2 and t = .1 is 

'^.05. It will be a bit larger for lower values of r. 

A part of this work ^2 was done in collaboration with 
Dr. Saurabh h. Eindani. 



Chapter 4 


ASYMMETRY PARiMETER AS A PUJCTION OP TRANSVERSE 
MOMENT IM IN THE PRODUCTION OP DI LEPTON 
WITH polarized BEAM AND TARGET 


4.1 Introduction 

In the previous chapter we investigated the angular 

distribution of one of the dimuon in dimuon restframe in 

13 

order to study the role of first order QCD effects in the 

dilepton production. This was necessitated because of the 

12 

fact that the simple Drell-Yan model does not explain all 
the details of dilepton data. In this chapter we study the 
production of massive lepton pair produced in longitudinally 
polarized hadronic collisions. We calculate the asjrmmetry 
parameter 


m 


+ - 


hr 


dm dq.^ 


X) - m^ u (H-,(+)H 

dm^dqj ^ 


m 


2 

2 d^o 


, 2 ., 2 ^"1 
dm dq^. 


2 

( Ht ( +) H,( +) p V'2) + m^ ( Ht ( +) H ^( -h 


^ 2-2 
dm aqg, 

X 100 


( 4 . 1 ) 


as a function of rrj^ = study of Aj^ as a function of 

will help us in two ways. At high values of q^, (greater than 


the average intrinsic transverse momentum of partons < GeV) , 

the asymmetry Aj^ is solely due to QCD processes because at high 
values of q^ the Drell-Yan contribution to the dilepton production 
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is almost zero. Conseq.uerLtly the variation of A-r with 
computed via first order QCI) processes can be compared with 
experimental data without any ambiguity and it can provide 
a good check on the performance of the QCD effects in dilepton 
production. Secondly, the variation of with r^ being 
dependent on the model chosen for the distribution of polarized 
partons inside polarized hadrons, will serve as a check on the 
different models proposed for spin-dependent parton distribu- 
tion functions inside polarized hadrons. 

The' first order QCD subprocesses giving rise to non- 
zero dilepton transverse momentum are; (i) qq annihilation 
giving rise to a gluon and a virtual photon eventually decaying 
into lepton pair, fig. 4.1(a), and (ii) qg Compton type inter- 
action giving rise to a virtual photon which decays into a 
lepton pair, fig. 4.1(b). ¥e calculate the assrmmetry parameter 
At for four different models of spin dependent parton distribu- 
tion function in polarized nucleon for p-p collision and pp 
collision. The Drell-Yan subprocess does not contribute 
because we restrict ourselves to high transverse momentum of 
dilepton. 

A similar calculation for polarized pp collision has 

26 


been done by Keisho Hideka. 
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4. II DeriTation of the QCD subprocess Cross-section and the 
Exipression for the Asymmetry Parameter A j_ 

¥e consider the following process: 

H^(P^) + H2(P2) |J-"(k3_) + anything (4.2) 

and H2 denote colliding hadrons with momenta Pj_ and P2 
respectively, and k2 are the momenta of the muon pair. We 
will work in the hadron centre of mass frame. The QCD sub- 
processes producing lepton pair are; 

(i) iCPjS) +q(k,s') ^.g(p’,e') + ii'’'(k^) + \L{k^) (4.3, 

where q, q are the annihilating quark and antiquark with 
momenta, p, k and spin s, a ' respectively and g is the gluon 
produced with momentum p ' and polarization e'. This process 
is depicted in fig. 4.1(a). 


(ii) q(p,s) + g(k,e) q(p’,s') + |J.'^(k^) + p (k2) 

q(p,s) + g(k,£) -»■ q(p’,s') + p'*'(k^) + p (k2) (4.4) 

Pig. 4.1(b) depicts this process in which a quark or an anti- 
quark of momentum p and spin s interact with a gluon of momentum 
k and polarization e giving rise to a quark or antiquark 
of momentum p' and spin s' and a pair of leptons. 

We define 


2 

and s, m > j s. 




= k^ + k2 
u, q^ and r^ are defined in (3.17). 


(4.5) 


^q^ f^^ y(x,+) denote the probabilities of finding 
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a quark q of helicity + with fraction x of the parent hadron's 
moment\im inside a hadron H of helicity + and - respectively 

%(_)*' ai*® respective probabilities 

for gluons . 

44{+)(=') = f^(+)(x,+) - 4(^)(x,-) 

f^(x) = 

=gH(4.)(^>+) Sh(+)(^»'”)=%( ( x,-) 

( 4 . 6 ) 

Parity invariance demands that 

gjj(^)(x,+) = gH(_)(x,-) gH(+)(x,-) = gjj(_)(x,+) 

and thus 

= Af|(x) 

Agjj(+)(x) = - Agjj^_j(x) = Agg(x) ( 4 . 7 ) 

The process q(p,s) + q(k,s*) ■♦•g(p'e') + p'*’(k^) + 

The amplitude for this process according to the diagram 
4.l( a) is 
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M = g [7(k,s') [t** y'' + Y**] u(p.s)] 


+ X ®' -,1 


(P')[x"'|-x] [uOcp) Y„ T^Ckg)] 
Q, 


( 4 . 8 ) 


After mod sq.uaring M and then introducing the spin projec- 


tion operators 


1 +J^ 1 + Y 

2 


2 — for the initial spins, we 
sum over all spins, polarization and colours and average 


over the initial colours and get 
iMi^ = (f) i-g^e^e^/q^) 1 


( 4 . 9 ) 


where 




( 4 . 10 ) 




•pv 


Tr 




-q' " 2p-p’ ’ " ' 2p’ - k ’ ^ ' 2 


ly i ri :. iL=V:. y ^' )1 

2p-p' ^ ^ ^ 2p'*k 2 ^ 


q‘ 

( 4 . 11 ) 


By straightfoi’ward calculation we arrive at the differential 
cross-section 

(q(p^)q(p ).».^^^i"x) = |(l-n^P2^^ 3^-^) ^ (i - - ±g - - 2) 

dm dt ^ s^ tu 

( 4 . 12 ) 

^2 ^^0 heli cities of the quark and the antiquark 

respectively, other symbols have already been defined in 
Chapter 3. Using the identity q^ = u t/s, we deduce 

-,2'^ . ic (l-Pi 

dm^dq^ ^ sV^(s-m2)2-4sq| 

^2 4 

- 2 ) 


( 4 . 13 ) 
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It is to be observed that the quark and the antiquark 
of the Same helicity decouple. [The helicities are defined 
in reference to the momentum directions of the individual 
particles. ] 


The process q( p, s) +g( k, e) q(p ' , s ' ) +p.'’’( + ^”(^2) 


The amplitude for this process depicted in fig. 
4.1(b) is 


M = 


g^e^e^ [S(p',s') [y^ ^ y'" + Y'^^^ Y*^] ^(p.s)] 


£v^k) -[x^ Xl • ^ [u(k3_) Y^v(k2)] (4.14) 

¥e mod square M and introduce the spin projection operator 
1 + 

for initial quark spin and then sum over initial and 

final spins and colours and average over the initial 
colours. ¥e get 


2 4 2 
• 1 , S e e. 


ImI^ = (h { — 




where L is the same as given by (4. 10) and 


( 4 . 15 ) 


= £;(k) £„(k) Tr ^ Y*") 




‘q'"' 2p*k ’■ ’■ 2p’-k 



2p‘k ' ' 2p’-k 


2 


( 4 . 16 ) 
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A straightforwaxd calculation gives 

2 2 

1^2 — ^ (q(T})g(A) [_ ia_+s.,.+ 2m _t l 
dm^du 93 ^ s5 


A p ^ p 

, . rU-2m s-2m ■, n 

+ nA {— ^ } ] 

s u 


(4.17) 


where r] and A axe the heli cities of the quark and gluon 
respectively. 

Using the identity = u ^ /s we get 
j^2 _ d^ _^ (q(q)g(A ) -♦• p'*’}a~X) = -2^ 


9sV( s-m^)^-4sq| 


[ 2 i 2 s! + . 2mfgzsh)+ {. ti±2s!l + 




s 




(4.18) 


Expression for the asymmetry parameter Aj^ 


To write an expression for the asymmetry parameter Aj^» 
we need the convoluted cross-section m — 2 * — 2 — (^^^(iij) 


dm dq| 


H 2 (h 2 )'^ p p"!) due to quark- antiquark annihilation suhprocess 

o .2 Uomp. + _ , 

(fig. 4.1(a)) and m — p — 3 (H 2 ^(h-^)H 2 (h 2 )-^ p p X) due to 

dm dq^ 

quark— gluon compton process (fig. 4.1(h)). Using equations 
(4.13) and (4.18) we get the following expressions for these 


cross-sections . 
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o j2_Annih l6a a 1 1 

in2 do — _ (h^)- ^Vx) = —27“^ ^ ^ ^^^2 

dm'^dq^ o o 

2-:; 

0{(s-m^)^ - 4sq|,}- V ^ [m^ ^ (q(Tl 3 _)q(n 2 >^Vi ^~^)] 

q=u,d,s q]_jn 2 dq^. 






2 


2 2' 




(4.19) 


^2 (S|^(h^)H2(h2)..n'"li X) = ^ f iXg 

dm dq^ ° ° 

e((S-m2)2 - 4Sa2j I I t,2 4!S (q(,)g(A). ^Vx)] 

q=u,d 5 S r\,h=± dm dip 




(Xi.ri)} 8 h fh )(=‘2'''> * 


^2^ 2 


+ 4 , n(x ,q)} g (xp,A)] 


2' 2' 
( 4 . 20 ) 


.lVx.p. 


where h^ ard h 2 are 
respectively. 


the helicities of the hadrons and H 2 


Using equation (4.19), (4-20), (4.6) and (4.7) we arrive at 
the following expression for the asymmetry parameter Aj^. 


r, ^ 1x100 ( 4 . 21 ) 

Ap - [ ^^nih_|_ 0^0 mp i^Annih ^ omp 
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^Annih 


-Comp 

}j 


Annih 


= m 


2 


m 


o a2^Annih 
Z 0. o 


^ ^ (H. (+)H (+)-► p ^l“x) - ^ - 0-^0 

dm^dqj ^ dm2 ^^2 


( H^( +)Hp(-)-^ p'‘’p x) = / dx^ J dx ,9 { ( s-m^) 2_4sq2} 


‘•'1 J 

O O “ 


^ 2 ^ 2 (q(+)q(+)‘^ P-'^P' x) -m"^ P 

dm dq^, dm^dqj 


(q(-t-)q(-)- p p-X)> [ I e2 U4 ^j(x^) 4- 

— u J Cl j S X 


(q-^-^ q)} ] 


"2' 
(4.22) 


.2 d^a^'^^P 


m- ( H. ( +) H-( +) p p X) - 

dm^dqj , ^ 

1 1 


m 


2 a2„Coni, 

2 2 
dm^dqj 


(Hq(+)H2(-)**> p^p”X) =1“ dx^ J dXpS { ( s-m^) 2_4sq^ 


2^ 2^ 

{m^ do (q(+)g(+) pVx) - 

dm dq^ dm dq^. 


(q(+)g(-)-^ p p X)} [ I \ I ^ Sg (+)(x^) ^ (+)( 2 ; 2 ) 


q=u,d, s 

- z, X 

+ Af^ (^ 2 )^ + ^ 

^2^^^ Hq H 2 


2 ' 

(4.23) 


„ .2^ Annih 

2 d-a— *“ (. (^,)H,(+).^.VX) + -"4-%-^ 


,2^Ainih 

... d O 

=z m 0 9 ^ ‘ ”^■^2 

dm'^dq^ 


dm dq^j 


{EA+)EJ-)-^ pVx) = /dxq J dX26{(s-m2)2 _ 4-sq2} 




{ m^ ^a-" 2 (q(+)q(+)-*‘ p'^p x;) + *^2^ 2 

dm^dq^ 

(q(+)q(--)^ pV~x)} [ I ^ 

q=u,d, s i ^ 


( A 7A\ 
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Comp 


P ,2 Comp 

m ^ ( H^( +) Hp( +) -*■ 


2 da 


2 Comp 


2 2 

dm'^dq^ 


(%{+)H2(-)- iiVx) = ! taj I (s-m2)2_4Sq2) 


{n,2 {(l(+)g(+)- nVx) + 

dm dq^, dm'^dq^ 


"■(<l(+)g(-)- 


nVx)}[ I e‘[gg (X3_){f| (Xj) + f| (Xj)} 
q — " 0 . ^ d 9 s 1 2 2 


X, X 

+ { ^ 2 }]] 


(4.25) 


4.III Spin Dependence of Parton Distribution Punctions 

We take the following four different models for the 
spin dependent parton distribution functions in longitudinally 
polarized nucleons. 

A. ¥e assume the sea and the gluon to be unpolarized and 

take the spin dependence of quark distribution f(x,+) based 
on a nonr elati vistic SU(6) wavefunction of the nucleon: 




where the second superscripts v and s denote valence and 
sea contributions respectively. Then 
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Afp(x) 

<3 

II 

= 1 

5 p 


Afp(x) 

= Af|(x) 

= - 5 


Afp(x) 

= 

= AfS(x) = if|(x) = 

Af^( x) = Af^ ( x) =0 

Agp(x) 

= Agp(x) 

= 0 . 

(4.27) 


We call it SU(6) model. 


B. Sehgal’s model : In this model we use the Sehgal's 

parametrization for Af'^(x) and Af‘^(x). Since for a nucleon 

P P 


<J > = <S > 1 + < S > -, + <L > 

z z quarks z gluon z 


(4.28) 


We take < L > = 0 . 

Zi 

•^^z^uark I [Af'^(x) 4A f^(x) +Af^(x) + Af‘^(x) + Af®(x) 

+ Af®(x)] (4.29) 

(^) = ;dx [Af'^Cx) + Af'^Cx) -Af'^(x) - a/(x)] (4.30) 

n.4.p 

=:/dx [Af'^(x) + 4f^(x) - AfS(x) - Af®(x)] • (4.31) 
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Thus 


<S > 




z q,uark 2 


G . 

+ 2(7^)= 

^17 


^ 0 ] + ^ /dx[A f®(x) +Af^(x) ] 


(4,52) 


Neglecting the strange quark contribution, one gets 

1 


where 


< S,>q^uark = ! - D) 


(7;^) = F + D 

'' Gy^n-*-p 


(4.33) 


(~) 

^Vh--ho 


= F - D 


Using the experimental values of D+F = 1.25±,01 and 

F/D = 0.5 +0.03, one gets < = 0.3+0.03. 

Hence 

<S > ^ = 0 . 2 + .05 

gluon - 

Assuming the sea to be unpolarized, we get 

0-3 = quark = I / +Af*(x)] 


(4.34) 


(4.35) 


The Bjorken sum rule is 

1.23 = I dx [Af'^(x) -Af^(x)] 


(4.36) 


o 


Assuming Af^x) and Af‘^(x) to be proportional to f (x) and 
respectively, the equations (4.35) and (4.36) give 
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aAx) = 0.456 

Af^(x) = -0.315 ( 4 , 37 ) 

¥s paX'ame'fcrize llie spin dependenl; gluon distribution as 

Ag(x) = px g(x) 

Such that (4.34) is satisfied. 

= j xp giz) ax = ;\i.x)”ax = p /2 

and hence p = 0.4 and we get 

Ag(x) = 0.4x g(x) 

g(x) = (an)(i.^)">/j, (4.38) 

We call this Seiigal’s model. 

C. Cahlitz and Kaur's model ; In this model, valence quarks 

interact with the sea and hence theix* spins are diluted. Let 
2 

sin 9 be the probability that a valence quark's spin changes 
in interaction with the sea. Suppose H(x) be the probability 
of Spin flip interaction between the valence and sea and 
H'(x) be the density of the sea relative to the valence. 
Suppose the number of particles having spin s before inter- 
action be 1, then the number of particles having the same spin 
will be 1 +1 If(x) H(x) and the number of particles having 
the opposite spin will be ^lT(x) H(x) and hence 


't 
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2 

sin 0 


N(x) H(x) 


(4.39) 


2S 

Gai’li'fcz sn.d. Kaun assuming ijiie sea "to be unpolaxizedj deduce 
the expression 


H(x) ll(x) = H^(l-x)^x"^ ( 4 . 40 ) 

The value = 0.052 is set hy the Bjorken sum rule, 

Eq. (4.36). A measure of the spin dilution induced by these 
interactions is given by 


Cos 20 = [H(x) N(x) + 1]"^=[1 + Hq( 1 -x )2 ( 4 . 41 ) 

The spin dependent structure functions axe simply given 
as the product of a function describing the asymmetries of 

29 

valence-quark spins as given by broken SU( 6 ) model of Close 
in the absence of the interaction, -with the spin dilution 
factor Cos 29. The structure functions for scattering off 
longitudinally polarized nucleons thus have the form 

2 g 5 (x) = 0 os 2 e [| fp{x) - I fp(x)] 

2g5;(x) = Cos2e [| fj(x) - I f®U)] (4.42) 

Since 

2s5(x) = I Af"(x) + I Af^(x) 

2g“(x) =|4fp(x) +iifp(x) (4.43) 

Eqs. ( 4 . 42 ) and (4.43) give 
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Afp(x) = Gos2e [fp'^(x) - I 

Afp(x) = Gosae [- ~ fp^^(x)] (4.44) 

The < Jg. >sum rule that the spin carried by the nucleon is | 
indicates that 11. of proton helicity is due to gluons. 

We parametrize the spin dependence of gluon distribution as 
before j i.e. 

Ag(x) = cx g(x) = c(^) (l-x)“ (4.45) 

and fix c such that 11,6^ of the proton's helicity is carried 
by gluons. ¥e get 

Ag(x) = 0.058 (m+l)(l-x)“ (4.46) 

D. Siver's model ; This is a model of spin dependent 

distribution functions in nucleons given by Babcock, Monsey 
30 

and Sivers. In this model they assume the sea to be 

polarized and parametrize the sea distribution, based on 

perturbation theory diagram in QCI) and the generation of the 
31 

sea? as 

f®^J^(x,+) = c fp®®'(x) [2 + (1-x)^] 

fp(?)(^’-) = ° [1 + 2 (l-x)2] (4.47) 

where c is to be fixed by the amount of momentum carried by 


by the sea. Then 
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= cf®(x) x(2-x) (4.48) 

They parametrize the gluon distribution, on similar considera- 
tion, as 

gp(4.)^x^+) = k g(x) [2 + (l-x) 2 ] 

gp(+)(xs-) - k g(x) [1 + 2(l-x)^] (4.49) 

Hence 

Ag (x) = kg(x) x(2-x) (4.50) 

where k is to be fixed by the fact that 50^ of the nucleon's 

momentum is carried by the gluon. Por g(x) of the form 

(l-x)^/x, k comes to be 

i, _ 

12rm+2) 

g(x) = (l-x)°^/x (4-51) 

Af^^(x) and Afp’^(x) are fixed by Bjorken sum rule (4.36) and 

the <J > sum rule (4.28). If one assumes the spin averaged 

s 

parton distribution as given in ecLUation (3.38), we get 

Af®(x) = a/(x) = Af'^(x) = 0.028(1-x)'^(2-x) 

P P P 

Af;j^(x)= 0.456 f^^(x) 

P P 

Af^^(x)= -0.315 

ir Jr 


(4.52) 
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|.IV Results ao-d Discussion 

¥e compute the as^mimetry pajrameter At as a function of 

gcalsd. tpansveise momentum = q.^/Y's fop the longitudinally 

polarized ppoton-antippoton collision and ppoton-ppoton collision 

foP the ■various models of spin dependent papton distpihution 

functions given in section 4. Ill fop fa = 27-4 Gev. Pop the 

spin avePaged paPton distpibution we take the ones of PeiePles, 

23 

Epueman and Wang given in equation (3.38) and the spin 
averaged gluon distribution is taken of the form k(l-x)“/x 
where k is to be fixed by the fact that 50 ^ of nucleon's 
momentum is carried by the gluons. ¥e do our calculations for 
m = 5 and 8 . 

Pop p~p collision, the gluon contribution is understandably 
very much less than the valence quark- anti quark contribution. 

We make the following remarks regardiag p-p collision. 

(i) Pop all spin dependent parton distributiong is negative 
for all values of Pjj . 

(ii) Pop the distribution obtained from SU( 6 ), i Aj^j increases 

as Prj, increases, the value of Aj^ ranging between -30^ 
to - 44 ^. 1 ' he variation of Aj^ with r^ for this model is 

sho-wn in- fig. 4 . 2 (a) and 4 . 2 (b), the gluon distributions 

g(x) being 3.0 (l-x)Vx and 4 . 5 (l-:ic)®A 

32 

The Drell-Ian model (with r^ = 0) also gives - 44 ^ when 

SU( 6 ) is used. This limit can be understood simply 

. ^ ■ ■ 2 \ 

because the u quarks dominate Aj^ and Af (x) = 5 


Thus 



Tr=j02 
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4 , 2 X ZJxv, X 

Q*X f (x) 




4 

9 ‘ 


(iii) Por Sehgal's model the parameter is more or less 
similar to SU(6) model but the values are less. The 
results are shomi in figs. 4.3(a) and 4.3(b). At 

Jj 

varies from -12^ to -20^ as r^ goes from .05 to .45* 

(iv) The Carlitz and Kaur model has a much larger value of 
Af^(x) for large x. This leads to a much larger 
asynmetry for large x as shown in figs. 4.4(a) and 4.4(b). 
Even though the model is not very reliable in this limit, 
it will be of interest to test where the asymmetry is 

as large as 80% for r^, = 0.45. This would clearly 
distinguish it from the other models. Aj^ varies from 
-30^ to -80%, as r,^, goes from .05 to .45. 

(v) Siver’s distribution, when the sea is also polarized, 
does not lead to any significant difference from the 
Sehgal's distribution. This is understandable because 
the contribution due to the sea is insignificant. 


¥e now turn to the p-p case where the quark-gluon 
contribution is dominant over the quark- antiquark contribution. 
This happens because anti quark distribution in a proton is 
suppressed. ¥e give the following compients. 

(i) For SU(6) Aj^ is zero because neither the antiquarks 
nor the gluons are assumed to be polarized. 
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(ii) Por Sehgal's distribution the asymmetry Aj^ is shown in 
fig* 4.5(a) for g{x) = 3(l--x)^/x and in fig. 4.5(b) for 
g(x) = 4.5(l-x)®/x. This comes entirely from quark-gluon 
scattering, as the' anti quarks ane assumed to be unpolarized. 

The asymmetry has a strong dependence on the gluon distribution. 
Por g(x) = 3(l-x)^/x, Aj^ varies from -5fo to +10^ and for 

O 

g(x) = 4.5(l-x) /x, Aj^ varies from vO to The same 

remark holds for the Carlitz-Kaur distribution. The varia- 
tion of Aj^ with r^^ for Carlitz-Kaur distribution are 'shown 
in fig'. 4.6(a) for-g(x) = 3(l-x)^/x and in fig. 4.6(b) for 
g(x) = 4.5(l-x)®/x. 

(iii) The antiquarks are polarized in the Siver’s model and 

hence both quark-antiquark and the quark-gluon processes 
contribute to the asymmetry. Por g(x) = 3(l-x) /x, Aj^ as a 
function of r,^ is shown in fig. 4.7(a) and it is in the 
range of ±% to for g(x) = 4.5(l-x) /x the variation of 

with r^, is sho>jn in fig. 4.7(b) where Aj^ is in the range 
of -5% to 

¥e have calculated the asymmetry parameter A^in the 
scaling limit. However, we have also repeated our calcula- 
tions by treating the parton functions as a function of 
q (=m ) and find no appreciable change in Aj^. ¥e use the 
Gluck-Reya distributions for the q dependence which is 
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(92/0.012) 

® P ° lil{ 

^av,„ „2,rln(92/7xl0-''),“-''91-=^l'" (l^/O.OOlS) 
^ ° ln(qj/7xl0 0 

f®(x q^) (<1^/0.583) 

^ ° ln(q^/lxlO“ )’^ 


9:.dV/ 2\ 

fp (x,q ) 


fp{x,q^) 


g(x,q^) 


ln(qQ/4xlO ^) 


with q' 


3 Gev"^ 


(4.53) 


¥e conclude by observing that the striking variation of 
Ap with r^ is a typical QCD phenomenon and hence is worth 
measurement in order to study the role of QCD in dilepton 
production. It also serves as a good test for the correct 
spin dependent distribution functions of quarks and gluons in 
polarized nucleons. 

This work has been submitted to Phys. Rev. D for 


publication. 



Chapter 5 


EXTRAGTIUG I' HE QUAjRX AND GLUON DISTRIBUTION 
IN MESONS PROM DILEPTON DATA 


5 .1 Intro duct,ion 

The parton model idea and the underlying QCD theory 
seems to work well in the asymptotic region of hadronic 
collisions. Por explicit calculations one needs the quark 
and gluon distribution functions inside hadrons. The 
quark distribution in nucleons are usually extracted from 
deep inelastic eN^ eX and vN-»- jiX data. The gluon distribu- 
tion is parametrized in the light of some theoretical 

25 

guidance such as quark counting rules, but there is no 

direct experimental check on the gluon parametrization. 

The data on inclusive production of massive dilepton in 

hadronic collisions can also be used to extract the quark 

17 

and gluon distribution. The usual extraction based on 
simple Drell-Yan model where gluons do not play any role 
cannot be very reliable. Considering the dilepton produc- 
tion in MN -»■ where M ia m— , K— , meson and N 

is proton or neutron, and using first order QCD perturbative 
effects at high transverse momentum of the dilepton we 
deduce in section 5. II certain relations which can be used 
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as a direct check on the quark and gluon distribution in 
meson, which requires only the vakence quark distribution 
in nucleon (i.e. the knowledge of sea quarks either in 
mesons or nucleons and gluon distribution in nucleons is 
not at all required). The basic assumption made is that the 
sea is an SU( 2) singlet . 

5. II Relations Between Processes X 

¥e will be restricting ourselves to high value of 
dilepton transverse momentum and consequently the simple 
Drell-Yan contribution from fig. 2.1 drops out. The non- 
zero contribution comes from diagrams 4.1(a) and 4.1(b) 
and the differential cross-section is given by 


m 


d^a 


2 2 
dm^dqj 




dx^ 


0{(§-m^) 



^ e^ [B (^ 2 ) 

+ c {%|(x3_)[f|(x2) + f|(x2)] 




+ [ 


^1 


M 


■H 




(5.1) 


where B and C are the differential cross-sections for the 

subprocess qq-*-gY* ^ _ snd qg.^ q^y* ^ _ respectively and 

■^ix p p. 

are given below 
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B = 


16 a 


27lv^(s-iii2)2_4-3q2 § 


^2 4 

[ — V - 2] 


(5.2) 


G 


a^a. 


9sf(s~m )^-4sq| ^ 


[ 1 - +^ ' "^ + {1 - 2a^(s-m^) 


>^2 

s 


>] 

(5.3) 


ill other symbols have the same meanings as defined in 
earlier chapters. 


We mil use equation (5.1) and the internal symmetries 
to derive a hierarchy of sum rules. The most reliable one 
will be the one based on SU( 2) . We assume that the valence 
quark distribution fpy^(x) of nucleons (p,n) and the 
^'^'^( 2 ;) of the meson doublets ( 71 ;“) , (K'^,K°), 

(K ,K ) are related through iso spin rotation and charge 


conjugation. 






fp^(x) = 

-dv/ >. 

II 


(1 

Hj 

= d"^(x) 

P 






( 5 . 4 ) 

f^;(x) 

n 

tC 

II 

II 

= 

% 

^Zix) 

% 


11 

K 

= fi;(x) = 

z 

fr(x) = 

Z 

ui(x) 

( 5 . 5 ) 

^x) = 

z 

Ox) 

z 

= ff!(x) = 

z° 

f®^(x) = 

Z~ 

Sg( x) 



We assume the sea quark distribution f'^®(x) and the 
gluon distribution g(x) in mesons to be an SU(2) singlet, i.e. 
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S ^(x) = g ^(x) 

"rt Ti~ 


Sjx) 


"Tt 


^ - S o(x) = g-«(x)=g (x) 

ii TT X E“ 




= gg;(x) 


n 

= f'^!(x) 

1 

II 

= t^lix) 

71 

= u®(x). 





f^^Cx) = 

= f^!(x) . 

= sj(x) 




% 

Ti: 

n ^ 

K 

= 

= f^o^x) 

= f‘*o(=') 
X 

= f^(x) 
£0 

= <"(-) = 


K" 

= f'^^Cx) 
iT 

= U®(x) 


- 

-pSS/ N 

f (x) 

K 

= f®®(x) 

K~ 

= f^(x) 

K 

= ff! (x) 

XO 

= s|(x) 

(5.6) 


Further the gluon distribution g(x) in nucleons is assumed 
to be Sir( 2) singlet and the up, down and strange luark sea 
are separately the same in proton and neutron 




^p(^) \ix) d®(x) _ d®(x) Sp(x) = s^ix) (5.7) 

The superscripts .r and s refer to valence and sea quarks 
respectively. 

Now consider the following six sets of reactions: 

1. Tc'^p 

^~P H- 
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2 - 

K“p ^ 

3 . Tt'n ^ 

4 . Z'n |a'^|a“l^ 
Z~n ^ \x'^ii~Y^ 

5* 2°P ^ tx+ia-Z^ 

K°p - 

6. K°n ^ txVz3 

2 °n - ^x"■^x-Z 4 


Using equation (5.1) 
sions . 


we can write the following t» exprea 


-dfaJJSi „2 d^cdU) 1 1 

■ dm^'/ "^1 / ■’^2ene-m2)2-4S,|] 






- - f’*'^(i )} +{f4’^(x 1 + fis, , , 


+ f^®(x 3 _) - - f 3 ®(x )}] 


1 ^ M 
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+0 [gijCzg) + f|®(xj^) 

■ - f«®(x)} 

M £ 

+ f%(-l) - gjj(Xi) ) {f|''(x2) + fWu^) + f|-(x2) 


+ f^^CXj)}] ] 

If we invoke equations (5.5) and (5-6) 


(5.8a) 


fj^ (x) _ f^^(x), f|^(x) = f3'^(x), 

ixi ri 

f3®(x) = f5®(x) 

M M 


(x) = f|®(x) = 


fjj (x) - f^®(x) 5 gj^(x) = g_(x), f|"'^(x) = 0, 


we gei: 


j,2 dfolMl , ji2 , r/* 2> 2 . ? 

= i i "^2® [(=-" - 4 s<i2] 

I el B [fl'^Cx,) - f5-(x.)] (5.8) 


q — u 5 d o s 


Biirdlarly j 


d^cr(]yin) ^2 d^a(Mn) p- } o o o 

^ Z2L2 = / J dx_9{(S-m2)2 _ 4-^2^ 


2. 2 


dm dq 


.2., 2 


dm dqj 0 


1 P 2 


4sq„} 




+ f|®(x^)} + {f|^(x2) + ff®(x2) - - ^P®(X2)> 


^f|^(xi) + f^®(Xp)} ] 


4- 
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0 [ {gii(^2) - Sj,(x2) ) { + fW(x^) + 4’'(x^)+ 4®(x^)) 

+ <0^2^ + fn®(^2> ^ + C(^2> " 


Using equations (5.5), (5.6) and (5-7), i.e. 


f^®(x) = fj®(x), f|®(z) = f^®(x), g^(x) = g^(x), we have 


P 






m2 -l^g(Mp) ^ ^ r 7- „2^ 2 

dm^dq^ dm^dq,^ 


dm dq„ 


dx^ J dx^9 [ ( s-m^) ^-4sqj3 


I e^ [B{f^^(x^) - f55^(x^)} {f37(xT) + f|®(xT)} 
q=u,d, s 


■q '- n '""2' *p '■"‘M ' "M ^ 1' 


■tc gi^(xj^) - *p'"(^2>’] ‘5.9) 

With the help of equation (5.8) and (5.9) we get the following 
relations . 


[„2 alg l K^ . . fniQi] / [x,2 4f|iK^ . a^q(K°n). ^ 

dm dq^ dm^dq^ dm'^dq^ dm'^dq^ 


_ rm2 .m2 

- Lm 0 0 dm^dqj 


dm^dq^ 


dm dq|p 


[„2 Ua( 4p ) . „2 l^°(gU) ]/|;„2 iUlIji 

dm'^dq^ 


2 2 
dm dq| 


2 2 
dm dq^ 


. „2 , d°(g°p) ] = 4 
dm dq^ 

(5.10) 

„2 U a L Efa l] = 1 

dm'^dq^ 

( 5 . 11 ) 


/ 
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2 2 
dm dq^ 


m2 

2 2 
dm'^dq^ 


/ dx-, J dXpS [(s-m2)2 




B- I [4p(x2) - 4uX(l,)] • u;^‘(x^) 




71 1 


(5.12) 


2 d^o(K%) 

P P 
dm^dq^ 


2 d^a{ 


2 2 ” •' 
dm dq^, o 


I dx. j dxpe [( s-m'^)2 _ 


4sq^] 


B[- I ^^^(^2^ 


(5.13) 


2 d a(m:^n) 
1 2 2 
dm dq^ 


2 d a{% n, 
1 2 2 
dm dq^ 


/ dx-, dx„e [(s-m^)^ 


4sq^] 


i [Up(=^2) - 4d^(=^2>l • 


(5.14) 


2 a^a(KV) 

iU P P 

dm dq^ 


-2 d2g(K~n: 

2 2 
dm dq^ 


= J dx-, J dXpG [(s-m2)2 _ 


^ P-i 
4sqm]' 


B[-fap(x2) u2(Xp)] 


(5.15) 


1 rm2 d^c( -n: p) 
c; L 2 2 

dm dq^ 




m2 ^.^1] + 1 [n,2 


2 2 
dm^dq^ 


2 2 
dm^dq^ 


.m2 ^2L^] 

dm'^dq^ 


= * / dx- / dx 8 [(s-m^)^ - 4aq2][u^(x ) . d^Cx,)] 


[C g^(x^) + B u®(xp] 


Ti;' 1' 


(5.16) 


1 r„.2 afpl g ^). _ „2 ^a (K°nl j ,-^2 dfa 


^ [m‘ 

^ dm^dq^ 


2.. 2 


dm dq 


_ m2 i-i^l 

2. 2 J 


dm dqj dm dq^ 


3- J dx^ J dx2e [{s-m^)^-Asq^][v7{xJ-d^{x )][C g^(x.) 






+ B u® ( X. ) 1 


.rv'l 
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5. Ill Di 


sous Sion 


These relations, especially (5-10) and (5.11), are not 
^rely a consequence of isospin transformation. One needs 
^oth isospin transformation and the explicit QCli expression 
(5.1) for these relations to hold, iin experimental check on 
such a relation thus is not a check on isospin transformation 
alone hut also on the QCD expression. 

The right hand side of the expressions (5.16) and 
(5.17) depend only on the sea and gluon distribution in meson 
and the valence quark distribution of the nucleon. The contri- 
bution due to the gluon is expected to be far more than the 

ibution due to the sea. Thus it could be a very reliable 
check on the various parametrisation of the gluon distribution 
in .uesons. Similarly either of the equations (5.12) and (5.14) 
can be used as a test for the parametrisation of the valence 
ribution u^(x), and either of the equation (5.15) and 
(5.15) leads to the determination of the nonstrange valence- 
quark distribution uj(x) in kaons. ThJ.s determination does 

p.nd on the parametrization of the sea and gluon distri- 
butions in either of the incident particles. The above method 
will be superior to the present determination^^ using only 
the Drell-Yan expression and neglecting the sea. 

An additional assumption about the sea and the gluon 
distribution that they are SU( 3) singlets, l.e. gjx) = gjj.(x). 



90 


_ Ug(x) = S®(x) = s®(x) and ^p(z) = d^Cx) = s®(x) gives 
lie following relations: 


dm^dqj 


2 d‘^a(x~-D) d^o(K p) 1 , r/* 2 n2 . 2-, 

‘ "J.J - ““ ~ :L2..2 = 9 I <^1 ‘ a=^2e [(s-a‘^)‘^ - 4sq^] 


T 


[B{[u^(xp) - ■^£(2:3_)][4 u^(x 2) + Su^Cx^)] + [ ] 
u^Cxp)} + C gp(x2){ 5[u^(xp) - Ug(xp)]+[uJ(x3_)-sJ(xp)]} ] 


(5.18) 


2 d^o(Tt"n) „2 d^a(K“n) 1 ?- ^ _.2^ 2 .^.2n 

1, ^L. _ m ^ ^ = 9 i ^^1 J ^^2® ((s-m ) -4sqj] 


2. 2 


dm dq^i 


dm dq^. 


[B(tu:^(x^) - uJ(x^)]-[4d^(x2)+5u®(x2)]+[v,J(Z;^)-sJ{Xj^)] 


u 


(x^)} + C gpCxp) { 5 [u3’(x^ )-uJ(x-, )]+[uJ(xt )-sJ(x-, )]> ] 


p.^.2 






(5.19) 


> d2cr(K°p) ^2 d^odc'^n) _.2 d^o(Ti^p) ^2 d^a(Tc^n) 

22~ 22“°^ 22 
dm dqj dm'^dqj dm^dqj 

+ 1 / dxp J dXpe [(s-m^) ^-4sq^] B[ u^( x^ ) -uj( x-, ) ] 


I 2 _/' TX - O . 


, 2 _,_+_ 


,2 


2 2 
dm dq^ 






(5.20) 


one further assumes SU( 5) then u^(x) = u^(x) = s^{x) and one 
ts the following relations from equations (5.18), ( 5 . 19)5 
.20) and (5.12), (5.13), (5.14) and (5.15): 



m 


m 


2 d a(Z^p) 

2 2 
dm^dqj 


m 


l^g(K+t) 

2 2 
dm'^dqjj 





91 

d^o(E“p) 

= 

2 d^a(TC~p) 

2 2 

(5.21) 

dm^dqj 


dm dq^ 


d^a(K"n) 

.2,2 



in2 d^g( m~n) 

2 2 

(5.22) 

dm'" dq^ 


dm'^dq^, 


_ 2 a^g(Z°n) 

2 2 
dm^dq^ 

_ m2 ^^2l£pl . m2 
dm^dq,^ 

(5.23) 

dm'^dq'^ 

= 16 m^ 

.2 

I^lsl 4m2 d^g(Tt~n) 

..2 .„2..2 

+ 4m2i!4£|i 

/V ^ 


dm dqj 


_ ni2 (5.24) 

dm'^dq^, 

l2_/Tr+„^ ^ j2_/'„— j2_/ +. 


- 2 dfoiKZ^ ^ 2 df^lnj, _ , 2 ^ 

)m 2_2 , 2,_2 . 2^„2 ^ .„2,„2 


dm dq^-, 


dm^dqj 


_2 d‘^g(7i p) 

dm'dqj dm^dq^^ 


. m2 d^g(%~ nl 
XU 2 2 

dm^dq^ 


(5.25) 


5 also define the following parameters 


= [m2 

dm dq^j, 


[m2 a^°(^P) 

dm'^dq^ 


m2 mlqi£ |l] /[m2 ^g i £ |l 

dm'^dq^ dm'^dq^ 


m2 a"°irP l]/[m2 

dm'^dq^ dm'^dq^ 


- m2 


2 2 
dm^dqj 


(5.26) 


. m2 


2 2 
dm^dq^ 


rrr,2 d^a(K'^n) 

— p p 

dm dq^ 


, m2 4!2^y[m2 

dm'^dq^ dm^^dq^ 


m2 a"°U°Pi .] 

dm^dq| 


(5.27) 
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9fr. 

For large r^ = the integrsoids in (5.12), (5.13), (5.14) 
and (5.15) contribute mainly near x^, ^ 1. I'hen in this 

limit one gets, 



1-4 u^(x)/d^(x) 
u (z)/d (x) - 4 

Jr J:/ 


Ej. Uj,(x)/dj,(x) 

u 

ibr SU(6) nucleon wave function 1^2 ~ ~ ^ 

large rgii whereas for nucleon quark distribution as given by 

'irueman, Peierls and Wang,^^ 1^2 = ® , Sj = -4- 

X-*- 1 

Hence and ^2 serve as good indicators for the tail distribu- 
tions of u (x) and d (x) . 

P P 

This work has been submitted to Phys. Ilev. D for 


publication. 



Ciiapter 6 


COUGLUaiOU 


12 

The Drell-Yaxi mechanism, in-volving on-shell q.uark- 

antiquark annihilation into a massive virtual photon, for 

lepton pair production in hadronic collisions does not 

explain all the details of the dilepton data.^"^ One needs, 

in addition to the Drell-Yan mechanism, QCD effects such as 

Q.1 ^ SY* + _s explain the dimuon data parti- 

-► P p” -►P P~ 

cularly the transverse momentum and the absolute normaliza- 
tion of the dimuon cross-section. The contribution due to 

13 

the first order QCD effects is found to be of the order 
12 

of Drell-Yan contribution. 

¥e have studied the angular distribution of the dimuon 

in their rest frame using first order QCD effects given by 

diagrams 2.3( a)&(b). The angular distribution is of the form 

1 + A Oos^e + B bin 26 Cos<!) + C Sin^O Cos 2 ^. A, B and C 

2 

are functions of m , s and q^ the transverse momentum of the 
dimuon. ¥e restricted ourselves to large q,j, in order to 
avoid mass singularity and infrared divergence and to 
justify neglecting primordial transverse momentum of the 
partons. ¥e have computed A and B (C = ^(1-A)) in the 
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Grottfried- Jackson frame, i.e, in the dilepton restframe, 
angles are measured taking the beam as the z axis and the 
X axis being in the plane of the two hadron momenta for pp, 
pp, m'^p, Tc'p collisions at fs = 27.4 G-ev. For low ralue of 
- 9 r|i/Vs < 0.05, the angular distribution is of the form 
1 + Cos 0, identical to the Drell-Yan prediction, but for 
higher values of r,^ > 0 . 05 , there is a strong, dependence on 
both 0 and <|) , which shows a striking departure from the 
hrell-Yan prediction. Hence we suggest the study of the 
dimuon angular distribution at large transverse momentum in 
order to investiga,te the role of QCD in the dimuon production. 

Purther we have studied the asymmetry associated with 

dilepton production in longitudinally polarized proton-proton 

and proton-antiproton collisions, using four different models 

of spin dependent quark and gluon distribution functions. 

For p-p collision the asymmetry is quite large (greater than 

15^2) and shows striking variation with r^, = q,^/Vs and is 

strongly dependent on tne model of spin-dependent structure 

functions. For p-p collisions the asymmetry is not large 

(a, 5 ^to 1%) and is also not strongly model dependent because 

in all the models, we use, the sea and the gluon are not 

strongly^ polar-ized. However, if one takes the model of 

34 

Bajpai and Ramachandran, which we have not considered in 
this thesis, where the gluon and the sea are polarized 
strongly but in opposite directions, the asymmetry is quite 



95 


35 

large (''*20^) for pp collisions also and changes sign if 
the number of flavours considered is taken to be greater 
than four. And thus the study of the asymmetry parameter A-r 
defined in equation (4.1) as a function of r^ will serve 
two purposes^ (i) it distinguishes the first order QCD 
subprocesses giving rise to dilepton pair where Aj^ will be 
nonzero from the quark-meson scattering processes giving 
rise to dimuon for which Aj^ will be zero and (ii) and it 
serves as a good test for the spin dependent quark and gluon 
distributions in polarized nucleons. 

finally we have derived relations between dilepton 
differential cross-sections in meson-nucleon collisions; 
n-H, K%, E;°N, (U = n,p) collisions, assuming the sea 

and the gluon distribution to be SU(2) singlet. These 
relations lead to tests of the parametrizations of (i) non- 
strange valence quark distribution in mesons and (ii) gluon 
distribution in n and K mesons. The knowledge of the sea 
quarks either in mesons or nucleons and gluon distribution 
in nucleons is not required. 

The effect of neutral vector boson has been 
56 

neglected because the mazimum value of dimuon mass m 
considered is approximately 16 Gev (for ][ s = 27.4 Gev and 


T = 0.4). 
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